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Abstract 

We elucidate the structural transitions in a helical off-lattice homopolymer induced by crowding 
agents, as function of the number of monomers (N) and volume fraction (cf> c ) of crowding particles. 
At (ft c = 0, the homopolymer undergoes transitions from a random coil to a helix, helical hairpin 
HH, and helix bundle HB structures depending on N, and temperature. Crowding induces chain 
compaction that can promote HH or HB formation depending on <f) c . Typically, the helical content 
decreases which is reflected in the decrease in the transition temperatures that depend on <f> c , N, 
and the size of the crowding particles. 
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ume fraction (<f) c ) of large macromolecules such as lipids, ribosome, and cytoskele- 
l| in the cell interior, which can be as large as 0.4 {2], affects all biological processes 



ranging from transcription to folding of RNA and proteins. Protein stability l3| and folding 
rates [5| of proteins are enhanced by an entropic stabilization mechanism (ESM) according 
to which the predominant contribution to the native state stabilization is due to an increase 
in the free energy of the unfolded states. Entropy decrease of the unfolded states results 
from the suppression of the number of allowed conformations of the polypeptide chains due 
to volume excluded by the crowding particles, while the native state is affected to a lesser 
extent. The ESM [3|, [5J is linked to crowding agent-induced depletion attraction 
between the monomers of the protein or RNA Crowding agents can also profoundly 
affect protein-protein interactions 10[ and amyloid formation 11] that is linked to a number 
of neurodegenerative diseases. 

We consider crowding effects on one the random coil RC to helix H transition. The 
interplay between the multitude of interactions between the crowding agents and proteins 
(Vcp) and the intra-protein forces (Vp) makes it difficult the structural changes that occur in 
a protein when <p c ^ 0. We consider the effect of spherical crowding agents on an off-lattice 
model of a homopolymer chain fl4|. which undergoes a coil to helix transition as temperature 
(T) is varied when C = 0. The major results, which were obtained for polymers with different 
N and <p c using molecular simulations, are: 

(a) The phase diagram is determined by a balance between the strength, 7, of the dihedral 
angle potential that is related to the local stiffness and the parameter 5, which specifies the 
strength of the hydrophobic attraction between the non-bonded beads, i.e. ones that are 
separated by three or more covalent bonds. As 5 is varied the homopolymer undergoes a 
series of structural transitions from a RC to H, helical hairpin (HH), and helix bundle 
(HB) at low temperatures, depending on N. For a fixed </> c , crowding particles whose radius 
(r c ) is commensurate with the the size of the monomer(r m ) [?], ~ 1), have the largest 
effect in stabilizing the collapsed structures. 

(b) The r c and </> c - dependent transition temperature, T$(r c ,(j) c ), from random coil (RC) 
to predominantly helical conformations (HH or HB) changes dramatically depending on r c 
and (j) c . For a fixed C , the most significant change in Ts(r c , <p c ) (compared to C = 0) occurs 
when r c is on the order of the size of the monomer. The values of Ts(r c , <j) c ) saturate, at all 
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N, when r c becomes large. At a fixed r c , Ts(r c , <f> c ) decreases as C increases. 

The conformations of a homopolymer chain, with N connected beads [ijj], The chain 
conformations are specified by the vector {r,} , % = 1,2,..N. The potential energy of 
the chain is a sum of bond-stretch potential, bond-angle potential, interactions associated 
with the (N — 3) dihedral angle degrees of freedom (Vb), and non-bonded potential (Vn) 
that determines the extent of tertiary interactions. The energy functions Vb and Vn are 
Vb = Eiiz {(1 + cos (<)), + f )) + (1 + cos 30;)} and V N = £\,. e 
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where e (=lkcal/mol) specifies the energy of interactions between non-bonded beads % and 
j separated by = \r, — rj\, 7 (=1 in this work) is the strength of the of the dihedral 
potential, 0, is the i th dihedral angle, and r m = 2 A is the size of a monomer. The potential 
y{— Vcc + Vcp), arising from interactions between the spherical crowding particles (Vcc) 
and with the monomers (Vcp) is V = £ 
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where is the 



j, r^- p is the distance between bead % and the 



distance between the crowding particles % anc. 
crowding particle j, and r c is size of the crowder. 

The simulations were performed using a modified in-house AMBER6 
was altered to incorporate Langevin dynamics in the low friction limit 



17l | package that 
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rate of conformational sampling 19j. In the presence of crowders the calculations were 
performed in the NVT ensemble. The homopolymer and the crowding particles are confined 
to a cubic box and periodic boundary conditions are used to minimize surface effects. The 
size of the simulation box, is determined by the condition that the box contain a minimum 
of 150 crowding particles. For small r c , the edge of the box is equal to the sum of the length 
of the fully extended helix and four times the average distance between the crowders at a 
specified <p c . The number of crowding particles ranges from m 150 for the most dilute system 
with the largest r c , to 1200 {<f> c = 0.2, r c = 2 A, N = 16), and 3000 {<f> c = 0.2, r c = 4 i, 
N = 64). 

The sampling; efficiency in the simulations are enhanced using the replica exchange method 
(REM) {5], [l^l , which ensures that the thermodynamic averages are fully converged. We used 
twenty replicas in the temperature range from T=100 to 400K in the REM simulations. The 
initial configurations are randomly chosen from high-temperature simulations, and subse- 
quently quenched to the desired temperatures. The integration time step is 10 _4 tl where 
T£ = (mr^/e)^ with m being the mass of a bead. At chosen time interval (= 40tl), config- 
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urations with neighboring temperatures are exchanged. The acceptance probability, which 
depends on the temperatures and the energies of the replicas, was in the range 0.2 — 0.3. 
In order to calculate averages we retained between (4,000 - 8,000) conformations for each 
replica. The results of the REM simulations were combined with independent data set gener- 



ated using the weighted histogram analysis method (WHAM) [20|] to obtain thermodynamic 
averages. The structural transitions in the hompolymer are characterized by using the spe- 
cific heat Cy and the radius of gyration R g . The extent of helical order was quantified using 
/j/(T, C ) = -— ^ J^T 3 (©(A0 — \(f>i — (f)f\)) where 9(x) is the Heavyside function, 0^ is 
the value of the i th dihedral angle in the energy minimized (T = 0) helical state, and A0 
is the tolerance in fa used to assign helical character to the i th dihedral angle. We chose 
A0 = 12.07° to ensure that Ts obtained using the criterion fn{Ts, C = 0) = 0.5 is consistent 
with the temperature at which Cy for N = 64 has a maximum. 

The structural transitions as a function of T and 5 for N = 16 with C = (Fig. la) show 
that the RC^H transition occurs at Ts ~ 292 K at 5 = 0. For low to moderate 5 values 
(6 < 0.5) the polymer exists either as a RC or a H. With 5 = 0.75, we find a transition to a 
HH that is accompanied by a drastic reduction in R g . At high 5 values, the energy cost to 
form a bend in HH is compensated by a number of favorable tertiary contacts that stabilize 
the HH. The chain compaction at high 5 results in structures that have high helical content 
as measured. For N = 32 and 64, besides H and HHs we find that helix bundles (HBs) can 
also form as 5 and T are changed. For some range of S the ordered structures coexist, while 
for other choices the probability distribution is peaked around only one unique structure. 
The transition temperature Ts at which the ordered structures form changes dramatically 
as 5 increases (Fig. lb). Only when 5 > 0.5, do we find significant dependence of T$(5, C ) 
on N (Fig. lb). 

From the temperature dependence of the thermally averaged R g with C = 0.2 and at 
various sizes of r c we find that smaller crowding agents (r c = 2 or 4 A) are most efficient in 
inducing chain compaction (Fig. 2a). For r c = 2 A the values of R g even at high temper- 
atures are considerably smaller than R^ - the radius of gyration of the energy-minimized 
helical structure (T = 0). A qualitative explanation follows from the Asakura-Oosawa (AO) 
theory [3], which predicts that the strength of the additional entropically-induced effective 
attraction between the beads increases as r r decreases and <b r increases. As a result, the 
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effective attraction 5r(0 c , r c ) ~ 5 + /(0 C , r c ), which increases with decreasing r c [21], is 
largest for small r c . In the helical homopolymer model a reduction in R g , at sufficiently 
large 5, is also accompanied by enhancement in helical order which explains the emergence 
of HH (Fig. 2b). Thus, crowding agents with ^ 1 are most efficient in inducing ordered 
structure formation at high T even if <5 is not large. 

The probability density, P(T, R g /R^), in Fig. 2b for <p c = 0.2, shows that the distribution 
function changes continuously as T decreases. At high temperatures, there is one broad 
peak that represents an ensemble of mostly random coil RC structures. As T decreases 
below Ts(r c ,(p c ) a sharp peak that corresponds to HH structures, with high /h(T, </> c ), 
emerges. Since there is only a continuous shift in P(T, Rg/R^), without a discernible region 
of bimodality, the transitions to structures with high Jh{T, <f) c ) are not "phase transitions". 
Rather, the energy landscape has multiple basins of attraction with varying helical content 
whose population can be altered by changing T, <p c , or 5. 

Transitions to higher order (three or more) HB structures occur for N = 64 with <p c = 0.2 
and r c = 4 A (Fig. 3). In this case, crowding- induced formation of HH and HB at low T is 
also accompanied by a dramatic reduction in R g (Fig. 3a). As temperature decreases, there 
is evidence for coexistence between long HH and HB (Fig. 3b). The formation of a large 
number of inter-helical contacts compensates for the energetic cost due to bend formation 
which results in the transition to the HB. 

The volume fraction C can be independently altered by changing either r c or the number 
density of the crowding agents. For a fixed N and r c , the values of Ts(r c ,<j) c ) decrease as 
C increases (Fig. 4a). The variations are larger for the smaller r c (Fig. 4a). The decrease 
in Ts(r c ,<f) c ) with increasing <j) c is a consequence of the enhancement in 5n(<f) c ,r c ) caused by 
the entropic depletion attraction. From the AO theory, it follows that the strength of the 
depletion attraction f((f> c , r c ) should increase as C increases. Thus, T s {r c cf) c ) should have the 
largest shift as <j) c increases, which is in accord with our simulations (Fig. 4a). 

The transition temperatures Ts(r c ,(p c ) ( obtained using /h-(Ts,0 c ) = 0.5) reports on 
the total helical content independent of whether the structure is a H, HH or HB. The 
changes in Ts(r c , 4> c ) for a fixed <p c and varying r c are shown in Figure 4A. We expect that 
as r c increases beyond the transition temperature Ts(r c ,(p c ) should approach the value 
expected for folding in narrow confined space formed by large crowding agents, and hence be 
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independent of r c . This is borne out by the simulations which show that Ts(<fr c , r c) is almost 
constant as r c > 30 A (Fig. 4b). For smaller values of r c Fig. 4b shows that T$ decreases 
sharply especially for N = 16. As in Fig. 4a, we find that the largest changes are obtained 
for r c = 2A. 

The stability of helical conformations is determined by interplay of the local stiffness 
and the specific attractive interactions between beads i and i + 3. In contrast, crowding 
agents, which enhance non-specific homogeneous attraction between the beads, induce chain 
compaction. Whether chain compaction is also accompanied by enhanced helical stability 
depends on the range and the strength of the AO attraction. At all values of r c , the transition 
temperatures Ts(r c , <p c ) decrease as r c decreases with the change being most dramatic for 
small r c (Fig. 4b). These results show that the helical stability decreases when C is non-zero 
even though the chain is compact. Thus, we conclude that the homogeneous AO attraction 
compromises the forces required to stabilize particular (H, HH, or HB) helical states. If the 
polymer backbone were stiff (7 > 1) then stretches of helical conformation on the scale of 

n 

the persistence length of the chain would persist especially at low temperatures [8( . For the 
flexible helical polymer there is a loss in helical stability, which is manifested by a decrease 
in Ts(r c , <f) c ) that is most pronounced when r c is small. 

We conclude with the following remarks. (1) Variation of the non-zero hydrophobicity 
parameter S leads to a variety of higher order structures at low temperatures, such as HH 
and HB, whose stabilities can be enhanced by macromolecular crowding. (2) The high 
temperature denatured structures are more compact at <p c ^ (see Figs. 2a and 3a) than 
when C = 0, which has implications for crowding-induced folding mechanisms of proteins. 
(3) The prediction that the helix stability changes, as C and r c are varied, can be validated 
using circular dichroism (CD) spectroscopy that detects the extent of helix formation in 
proteins 
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Figure Captions 

Figure 1. (a) Normalized radius of gyration R g /R^ as a function of temperature for 
several values of the short-range attraction parameter 5 for the N = 16 chain. Snapshots of 
typical conformations are also shown. The value of R^ = 10.1 A . (b) Transition temperature 
Tg as a function of 5 for iV = 16, 32 and 64. 

Figure 2. (a) Normalized radius of gyration R g /R^ as a function of T for different crow- 
der radii r c (indicated in the Figure) at C = 0.2 and N = 16. (b) Probability distribution 
functions of R g /R^ for different temperatures for <p c = 0.2 and r c = 2 A case in (a). 

Figure 3. Same as Fig. 2 except N=64. (a) Normalized radius of gyration R g /R^ (with 
R^ = 40.4 A ) at constant <p c = 0.2. (b) Probability distribution functions of R g j R g for the 
4> c = 0.2 and r c = 4 A case shown in (a). 

Figure 4. (a) Helix-coil transition temperature Ts(r c , <p c ) as a function of C at a constant 
r c for three chains with N = 16 (r c = 2A), N = 32 (r c = 4A), and N = 64 (r c = 4A). (b) 
Tg(r c ,(j) c ) as a function of the crowder radius r c at C = 0.2. Plots for chains of lengths 
iV = 16, 32 and 64 are shown. In both (a) and (b) 5 = 0.3. 
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